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SIGNAL GENERATION AND MIXING The different potential energy wells 10 and 12 represent 

ELECTRONICS FOR FREQUENCY-DOMAIN another type of energy in the diatomic molecule: electronic 

LIFETIME AND SPECTRAL energy. Different electronic energies occur when the mol- 

FLUOROMETRY ecule absorbs energy in such a fashion that it causes an 

5 electron to move to a higher energy configuration within the 

The U.S. Government has a paid-up license to certain molecule. The classical example of the energy absorption is 

technologies disclosed in this invention and the right in the 'changing the spring' that connects the two atoms. If we 

limited circumstances to require the patent owner to license "add" electronic energy to the potential energy well 10, the 

others on reasonable terms as provided for by the terms of energy raises the potential energy well 10 to the energy level 

the following contracts: NAS1-20426 and NAS1-20162 10 of potential energy well 12. 

awarded by NASA; DAAH01-91-R198 awarded by US With this background, we can use FIG. 1 to describe the 

Army Missile Command; DAAL06-92-C-0014 and process of fluorescence and fluorescence lifetime for a 

DAAD07-91-C-0127 awarded by US Army White Sands diatomic molecule. The example molecule is originally in a 

Missile Range; DAAA15-93-C-0034 and DAAM01-94-C- state depicted by the lower potential energy well 10. The 

0033 awarded by US Army Chemical and Biological 15 molecule absorbs a photon of the correct amount of energy 

Defense Agency; and F41624-95-C-6010 and F41624-97- that induces an electron to move to a higher electronic 

C-6029 awarded by US Air Force Armstrong Laboratory. energy state, represented by the upper potential energy well 

„ „ „ nn — m _ 12. The diatomic molecule also tends to undergo a change 

BACKGROUND OF THE INVENTION from , . lowef vibrational state > (the horizontal j^es 20 and 

1. Field of the Invention 20 22) in the lower potential well 10 to a higher vibrational state 

The present invention generally relates to fluorescence in the upper potential well 12. This initial excitation of the 

.U measurements. More specifically, the present invention molecule is shown by line 14. Through any one of numerous 

Issx relates to frequency domain measurements of the fluores- possible processes, the higher vibrational state 24 in the 

w cence lifetime and the fluorescence spectrum. u PP er wel1 deca Y s t0 a lower vibrational state 22 in the upper 

0 2. Description of the Related Art 25 ? eU (aS sh f wn by . L Une l *\ ^ P^cess of non-radiative 

1 s 2 ■ decay to a lower vibrational energy level in the upper well 
Siftf The process or fluorescence occurs when a substance, ' ... V~ . , 

. v . , , . A * « * occurs very rapidly. After the non-radiative vibrational 

y | such as a molecule, absorbs light at cue wave eng b or rekxati [ he ^Jon will want to return or revert back to 

£ energy), and then emite light a a longer wavelength or a ^ eQ ^ ^ showQ ^ u (k o(hef wof(J ^ 

* lower energy). A slight time delay occurs from when the , n . , °; . , , J . . . v , e m 

if* , / •_ . i_ .t_ i_ 30 electron wants to move back to where it was before the 

|J ■ substance absorbs light and when ttie substance re-emits M ^ ^ m ^ do[ ^ molecu , e emits 

|* hghl at the longer wavelength. Tto time delay is known as a ^ q£ ^ { J h a ^ of eQergy ^ ^ 

s| the Uuorescence litetime. excitation photon). This lowering of the energy state and the 

u FIG. 1 depicts the fluorescence process schematically for emission of a photon is the process of fluorescence. The 

; a diatomic molecule, which is a molecule that is comprised 35 process of emitting the photon has an associated time delay, 

f=- of two atoms. There are several "types" of energy that are which fe the fl uoreS cence lifetime. 

5~l internal to a diatomic molecule, one of which is vibrational . . ... . ,.~ 4 t c 

, . i / u i ~ - -~ In practice, we typically probe many different types of 

sx energy. The classic example of vibrational energy is imag- , r . * J \. ' r . . v . J S, ir> - 

:5= . . ^ . * j * »u. u • u molecules at once with the excitation light pulse. FIG. 2 

— ining the two atoms connected together by a spring where ... , 4 lL , • . i_ « „• i c 

e . n ♦ u i a r 7u i *u • e »u illustrates the case where a smgle short excitation pulse of 

the atoms oscillate back and . forth along the axis or. the 40 .. . . , . , . 7 e , , r , 

A . u u *u .u • hght 30 is absorbed by a sample of identical molecules all 
5 s; spring. As the two atoms approach each other, they expen- r « fl j 11 u- c 
I" 55 v & , # iU f *u • 1 at once. The fluorescence decay curve 32 resulting from a 
ence repulsion due to the proximity of their negatively . a ^ . c r . , 
, II, 1 j a fu . 11 <r typical fluorescence response of a sample of identical mol- 
charged electron clouds. As the atoms pull apart from one Jr . . 1 • * t_ .nr.uj 
another, they experience an attraction that one can imagine ecu es » exponential ui nature because not all of the .den- 
as the result of the attractive forces between the positively 45 teal molecU if em,tS lK Ayorescencc photon at precisely the 
charged nucleus of one atom for the negatively charged ^ «ponenUal decay foUows a mathematical 
. f , , r 4 . 4 1 » J ^. function so that we can calculate the fluorescence lifetime, 
electron cloud of the other atom, and vice versa. The m 9 

distance between the two atoms is known as the interauclear T * 
distance, which constantly changes as the two atoms oscil- 
late back and forth. As the interauclear distance decreases, 50 fy)*!^'-**" 
a rise in potential energy occurs due to the repulsive forces; 

as the interauclear distance increases, a rise in potential where ^ is the time the excitation pulse, I 0 is the initial 

energy again occurs due to the attractive forces. FIG. 1 fluorescence and I(t) is the observed fluorescence intensity 

illustrates this difference in potential energy by line 11 in the as a function of time. 

potential energy "well" 10 and by line 13 in the potential 55 l n the cases with multiple molecular types (where more 

energy "well" 12 where line 11 and 13 plot the potential than one lifetime decay is present within a sample of a target 

energy as a function of interauclear distance. mixture), the exponential decay seen in FIG. 2 would appear 

A group of horizontal lines 20, 22, and 24 appear in each as a sum of exponential functions. Prior art fluorometers 

of the potential energy wells 10 and 12. Each horizontal line typically use a type of time-correlated (or time-resolved) 

represents the individual vibrational energy states possible 60 system that count the emitted photons (from an excited 

for the diatomic molecule. Quantum mechanics requires that molecule) in order to measure the fluorescence lifetime, 

the vibrational frequency of the spring be within certain Other prior art systems add the ability to take the fluores- 

values. In other words, the spring may oscillate at frequency cence lifetime measurements in the frequency domain by 

'a', and the spring may oscillate at frequency 'b\ but it is a modulating the dynode of a photomultiplier tube, followed 

physical impossibility for the spring to oscillate at any 65 by a mixing and correlation procedure. These prior art 

frequency between V and V. In FIG. 1, any two adjacent systems are cumbersome, time consuming, and complicated 

horizontal lines would represent V and V. to operate. The present invention overcomes the limitations 
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of the prior art systems by utilizing a novel technique to FIG. 1 illustrates the fluorescence process for a diatomic 

measure the fluorescence lifetime and spectrum. Instead of molecule. 

irradiating the target sample with a single short pulse of light pjo 2 illustrates the absorption and emission of energy 

(photon counting), the present invention continuously irra- fof the fl uor escence process. 

diates the target sample with a light source whose amplitude 5 _ _ .„ , . . f 
modulation frequency is stepped with time. This technique ^ lllustrates lhe . chemometric extraction process of 
aUows us to use a chemometric analysis to automatically lhe P resent * v ™t">n for fluorescence spectrums and life- 
extract the lifetimes from the 'phase delay* and * intensity vs. times. 

frequency* characteristics of the emitted light. FIG. 4 is a block diagram of the present invention. 

SUMMARY OF THE INVENTION 10 FIG. 5 is a block diagram of a component of the present 

rpii . • . • c u i invention for generating an excitation signal. 

The present invention is a system for chemometric analy- & ° \ ° 

sis for the extraction of the individual fluorescence spectrum FIG - 6 is a block diagram of an embodiment of the present 

and fluorescence lifetime from a target mixture. The present invention for processing an emitted signal. 

invention comprises a processor with an apparatus for FIG. 7 is a block diagram of the present invention that 

generating an excitation signal to transmit at a target mixture * 5 illustrates the electrical pathways of the excitation signal and 

and an apparatus for detecting the emitted signal from the the emitted signal. 

target mixture. One embodiment of the present invention FIG. 8 is a block diagram of an embodiment of the present 

uses a processor that comprises a computer that extracts the invention for a chemometric analysis of the emitted signal to 

individual fluorescence spectrum and lifetime measurements determine the fluorescence spectrum and fluorescence life- 

from the frequency and wavelength data acquired from the 20 t j mes 

emitted signal. The present invention first determines the G nG 9 fc anothef block $ of an embodimen , of the 

and S matrices from the frequency and wavelength vectors. inventioQ for a cbemometric analysis of tbe emitted 

n A ^normalization of w S to G, occurs next due to the a j (0 determine tbe fluorescence spectrum and fluores- 

'=■. differences in units. Next, the present invention determines ° i :f - ( ;_„ 

*"» * . -i ™ ^ cence liletimes. 

IbJ the initial U, V, g and w J s approximations. The present & .„ .„ . , . . , , ... 

[ , 1 invention uses an iterative solution that first requires the p 10 lUustrales the P««nt mvenUon s ability to dlf- 

^ initiahzation of the decision variables including r,. The ferem ' ate m a tor S el mixture - 

•^J iterative solution compares the decision variables for con- DETAILED DESCRIPTION OF THE 

«§« verge nee to see if further approximation determinations of INVENTION 

ffl the U, V, g and w~*s are necessary. When the solution 30 

j\ converges, the present invention then determines the The present invention is a system for chemometric analy- 
reduced best fit error for the analysis of the individual sis for Ac extraction of the individual component fluores- 
5! fluorescence lifetime and the fluorescence spectrum. And cence spectra and fluorescence lifetimes from a target mix- 
Is* finally, the system of the present invention extracts indi- aire. Additionally, the present invention is a method and 
! . vidual component fluorescence lifetimes and fluorescence 35 apparatus for generating and mixing signals for frequency- 
C spectra from the emitted signal of the target mixture. domain lifetime and spectral fluorometry. This disclosure 
O The present invention additionally includes a method and describes numerous specific details that include specific 
iO apparatus for generating and mixing signals for frequency- processes, structures, and circuits in order to provide a 
n domain lifetime and spectral fluorometry. The present inven- lhorou S h understanding of the present invention. For 
^ tion comprises a signal generator that generates a driving/ 40 exarD P le ' lhe P resent invention describes the extraction of 
I"* reference signal that modulates the amplitude and/or the lhe fluorescence spectra and fluorescence life times from a 
frequency of the driving/reference signal over time and a tar ^ nnxnare However, the practice of the present inven- 
signal generator that generates the mixing signal that modu- tl0n mcludes f ^ extraction or analysis of other physical 
lates the amplitude and/or the frequency of the mixing signal Properties of a target mixture other than the previously 
over time. The driving/reference signal generator drives an <* described ones One skilled in the art will appreciate that one 
excitation signal generator that generates the excitation ™y P ra A c ^ e . lne V*** invention without these specific 
signal. The present invention then directs and transmits the Additionally, this disclosure does not describe some 
excitation signal at a target mixture, which absorbs the wel1 known processes such as Fourier transforms, 
energy from the excitation signal. The property of fluores- eigenequations, eigen functions, eigenvalues, transforms, or 
cence causes the target mixture to emit an emitted signal that 50 best fit analyses in detail in order not to obscure the present 
the present invention detects with a signal detector. The invention. 

present invention uses a mixer that mixes the emitted signal ^ P rescnt invention utilizes a novel system that uses a 
with the driving/reference signal to produce a processor chemometric analysis to determine the fluorescence life- 
reference signal. Another mixer mixes the emitted signal limes ^d the fluorescence spectra of a target mixture, 
with the mixing signal to produce a data signal. The present 55 Instead of irradiating the target sample with a single short 
invention then uses a processor to compare the processor P^^ of ^ (P hoton counting) as other prior art systems, 
reference signal with the data signal by analyzing the th e present invention continuously irradiates the target mix- 
differences in the phase and the differences in the amplitude ^e with a light source whose amplitude modulation fre- 
between the two signals. The processor then extracts the q ueDc y * sle PP ed with time. FIG. 3 illustrates the process of 
fluorescence lifetime and fluorescence spectrum of the emit- *o the present invention. The excitation light or signal 42, 
ted signal from the phase and amplitude information using which is the signal the present invention uses to continu- 
a chemometric analysis. ' ousl y irradiate the target mixture, is amplitude and phase 

modulated in a sinusoidal fashion. As previously discussed, 

DESCRIPTION OF THE DRAWINGS me emitted signal 40 is the photon emitted (emission 

To further aid in understanding the invention, the attached 65 intensity) from the molecule as a result of the fluorescence 

drawings help illustrate specific features of the invention and process. The emitted signal 40 will appear as a delayed 

the following is a brief description of the attached drawings: sinusoidal-intensity light of a longer wavelength. Embedded 
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within the emitted signal is the information necessary to 
extract the fluorescence lifetime and spectrum. However, 
one must use a chemometric analysis to extract the fluores- 
cence lifetimes from the 'phase delay <P' of the emitted light 
(compared to excitation signal) and the 'intensity vs. wave- 
length' (or spectrum) characteristics of the emitted light 
(compared to the excitation signal). 

FIG. 4 is a block diagram of the apparatus component of 
the present invention. A processor 50 couples to a data 
acquiring component 52 through a data input/output path 68. 
The preferred embodiment of the processor 50 comprises a 
computer. The data acquiring component comprises a gen- 
erator component 54 and a receiver component 56. Coupled 
to the generator component 54 Is an excitation signal gen- 
erator 58. The excitation signal generator produces an exci- 
tation signal (or light or photon) that is transmitted to a target 
mixture 66. The target mixture 66 absorbs the electronic 
energy of the excitation signal as the molecules of the target 
mixture move to a higher electronic energy state. During the 
fluorescence process, the molecules of the target mixture 
emit photons (an emitted signal 64) as the molecules move 
to a lower electronic energy state. A detector 60 detects the 
emitted signal 64 from the target mixture 66. Coupled to the 
detector is the receiver component 56, which processes the 
returned emitted signal 64 into a reference signal and a data 
signal for the processor 50. The processor 50 then uses a 
chemometric analysis to determine the individual fluores- 
cence spectrum and fluorescence lifetime of one or more 
fluorophores within the target mixture 66. 

FIG. 5 is a block diagram of the generator component 54 
of the present invention for generating the excitation signal. 
The present invention uses heterodyning techniques to pro- 
duce two sinusoidal RF signals, a driving/reference signal 
90 and a mixing signal 92. The present invention modulates 
the frequency of the two signals from 10 MHz to 200 MHz. 
One skilled in the art will appreciate that the present 
invention could vary the signals over a much larger fre- 
quency range. The preferred embodiment of the present 
invention generates the two signals with a frequency differ- 
ence of 10 kHz. Another embodiment of the present inven- 
tion generates the two signals using an adjustable offset 
frequency where the offset frequency is set through to the 
present invention's control software. 

Referring back to FIG. 5, the present invention generates 
the driving/reference signal with a driving/reference signal 
generator 71 and generates the mixing signal with a mixing 
signal generator 73. The driving/reference signal generator 
71 comprises an oscillator 74, a phase-locked-loop (PLL) 
70, a mixer 78, and a signal splitter 72. The mixing signal 
generator 73 comprises an oscillator 76, the phase-locked- 
loop (PLL) 70, a mixer 80, and the signal sputter 72. The 
present invention generates the two signals by mixing a 
signal from the PLL 70 with the signals from the two fixed 
oscillators 74 and 76. The preferred embodiment of the 
present invention operates oscillator 74 at 290 MHz and 
oscillator 76 at 290.01 MHz. One skilled in the art will 
appreciate that other operating frequencies for oscillators are 
possible. The preferred embodiment of the present invention 
can set the output of the PLL 70 to a frequency ranging from 
300 to 560 MHz, which the present invention controls 
through the control software of the present invention 
through the digital I/O datapath 68. 

Coupled to the driving/reference signal generator 71 is a 
filter 82 and an amp 86. The filter 82 filters out any unwanted 
components from the driving/reference signal, and the amp 
86 provides additional gain to the signal if necessary. 
Coupled to the mixing signal generator 73 is a filter 84 and 
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an amp 88. The filter 82 filters out any unwanted compo- 
nents from the mixing signal, and the amp 86 provides 
additional gain to the signal if necessary. The present 
invention uses the driving reference signal 90 as a phase and 

5 amplitude reference signal on the receiver component 56 (of 
FIG. 6). And, the present invention sends the mixing signal 
92 to the receiver component 56 where it is split. 

The present invention uses the driving/reference signal 90 
to directly drive and modulate an excitation signal generator 

10 58 (of FIG. 7) that may comprise a laser diode, a light 
emitting diode (LED), or a Deuterium lamp. The present 
invention can also use the modulation to directly modulate 
an electro-optic crystal for use in external modulation of gas 
laser sources or other continuous emitting excitation light 

15 sources. 

FIG. 6 is block diagram of the receiver component 56 of 
the present invention for processing the returned emitted 
signal. The present invention receives the modulated emitted 
signal 64 from the target mixture through a detector (60 of 

20 FIG. 7). The receiver component 56 receives the detected 
emitted signal 96 from the detector through an amplifier 97 
and receives the driving/reference signal 90 through an 
amplifier 99. The receiver component 56 additionally 
receives the mixing signal 92 through a splitter 98. The 

25 present invention then mixes the detected emitted signal 96 
with the mixing signal 92 with a mixer 102 and mixes the 
driving/reference signal 90 with the mixing signal 92 with a 
mixer 100. The resulting signals contain the sum and dif- 
ference frequencies of the mixer's (100 and 102) input 

30 signals. The present invention generates a data signal 114 at 
10 kHz by filtering the mixed detected emitted signal 
through a filter 106 with an amplifier 110 to eliminate the 
high frequency components. The present invention addition- 
ally generates a processor reference signal 112 at 10 kHz by 

35 filtering the mixed driving/reference signal through a filter 
104 with an amplifier 108 to eliminate the high frequency 
components. One skilled in the art after reading this disclo- 
sure will appreciate that one can use frequencies other than 
10 kHz to process signals. An analog to digital converter 118 

40 converts the data signal 114 into digital data for transmittal 
over the digital I/O path 68 to the processor (50 of FIG. 4), 
and an analog to digital converter 116 converts the processor 
reference signal 112 into digital data for transmittal over the 
digital I/O path 68 to the processor. The present invention 

45 then processes the two input signals to determine their 
relative amplitudes and phases because the amplitude and 
phase of the 10 kHz signals are directly proportional to the 
amplitude and phase of the high frequency modulation 
signals. As the modulation frequency is stepped from 10 

so MHz to 200 MHz, the processor reference signal 112 and the 
data signal 114 remain at 10 kHz. 

FIG. 7 is a block diagram of the present invention that 
illustrates the electrical pathways of the excitation signal and 
the emitted signal. The present invention generates the 

55 driving/reference signal with the driving/reference signal 
generator 71 and generates the mixing signal with a mixing 
signal generator 73. The present invention uses the driving/ 
reference signal 90 to directly drive and modulate an exci- 
tation signal generator 58 that may comprise a laser diode, 

60 a light emitting diode (LED), or a Deuterium lamp. The 
present invention can also use the modulation to directly 
modulate an electro-optic crystal for use in external modu- 
lation of gas laser sources or other continuous emitting 
excitation light sources. The excitation signal generator 58 

65 produces an excitation signal 62 that is transmitted to a 
target mixture 66. A photodetector 120 detects the excitation 
signal 62 on its way to the target and regenerates the 
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driving/reference signal 90 for retransmission to the receiv- line (linear least squares fit) to a graph of log(I(t)) vs time, 

ing component 56 of the present invention. The target we can determine the lifetime of a fluorophore. However, if 

mixture 66 absorbs the electronic energy of the excitation more than one fluorophore is present, the intensity variation 

signal as the molecules of the target mixture move to a with time depends on the concentrations, c„, and lifetimes, 

higher electronic energy state. During the fluorescence 5 x «» of the ^dividual fluorophores: 
process, the molecules of the target mixture emit the emitted 

signal 64 as the molecules move to a lower electronic energy /(/) = c n e' ( '-^ Wr " (4) 
state. A detector 60 detects the emitted signal 64 from the « 
target mixture 66. The detector 60 could comprise a photo- 
detector. Additionally, wavelength dispersion of the emitted 10 Nqw ^ d ^ mus( be fo by a non ^ Q ^ T 
signal can be accomplished by using a spectrophotometer, leas( procesSi which requ ires initial guesses for the 
monochromator or low-pass, high-pass and notch optical num ber of fluorophores, r, their concentrations, and their 
filters. From the detector 60, the present invention mixes the lifetimes. In general, convergence of the solution cannot be 
signals as previously described for FIG. 6 above. The guaranteed as it depends on the accuracy of the initial 
present invention then uses a chemometric analysis to deter- 15 guesses. 

mine the individual fluorescence spectrum and fluorescence Weber showed that decay constants and relative intensi- 

lifetime of one or more fluorophores within the target ties of an arbitrary number of independent components in a 

mixture 66. heterogeneous fluorescence emission could, in principle, be 

The present invention comprises several techniques or determined by measuring phase shift and relative modula- 

alternative embodiments to "target" the target mixture. One 20 tion of the total fluorescence at an equal number of frequen- 

embodiment of the present invention uses a target compart- cies. See Weber, G., Resolution of the Fluorescence Life- 

ment for holding the target mixture within a chamber. A times in a Heterogeneous System by Phase and Modulation 

container may contain the target mixture, or the mixture may Measurements, J. Phys. Chem., 85, 949-953 (1981). At each 

be free floating, in a fluid state, within the chamber. Another frequency, the phase and modulation measurements yield 

embodiment of the present invention uses an optical probe 25 the complex Fourier transform of the fluorescence impulse 

that combines the signal generator 58 and the detector 60 response, G+iS. The moments of a distribution of lifetimes 

into a single probe. A user of the present invention would are found as linear combinations of the G's and S's and 

then probe or place the optical probe into the target mixture. Prony's method is used to obtain lifetimes and concentra- 

And, another embodiment of the present invention may use tions of the components from the moments. Due to its 

a laser to target a free floating target mixture, such as the 30 sensitivity to noise, this method is really not practicable for 

emission from a car (at a distance). One skilled in the art will more than binary mixtures. 

appreciate that other embodiments of targeting the target Lackowicz and his coworkers developed an analysis of 

mixture are possible after reviewing the disclosure of the frequency domain phase and modulation data based on a 

present invention. weighted non -linear least square analysis. See Lackowicz, J. 

FIG. 8 is a block diagram of an embodiment of the present 35 R -> et al » Analysis of Fluorescence Decay Kinetics from 

invention for a chemometric analysis of the emitted signal to Variable-Frequency Phase Shift and Modulation Data, Bio- 

determine the fluorescence spectrum and fluorescence life- P h y s - J > 46 » 463-^77 (1984). By taking data at several 

times. Before discussing the block diagram of FIG. 8, we emission wavelengths, they were able to resolve binary 

first need to discuss the mathematical foundation used for mixtures whose components differed in lifetime by 30% and 

the chemometric analysis. Fluorescence occurs when a 40 ^raary mixtures the ratio of whose range of lifetimes 

molecule, excited by a photon to a higher electronic state, spanned a decade. Best fit is determined by minimizing a 

decays non-radiatively to a lower vibrational level of the rcduced chi-squared function. Poor initial guesses can still 

excited state, and then decays radiatively to a vibrational produce unpredictable results. 

level of its ground electronic state. The non-radiative decay Working with an Excitation-Emission-Frequency Array, 
is assumed to be 'fast'. Hie characterization of the radiative 45 McGown and her coworkers employed a principal compo- 
decay is the fluorescence lifetime, x. And as previously nents analysis followed by a non-linear least squares mini- 
stated, the wavelength of the fluorescent photon is greater mization on an over determined set of data to extract the 
(meaning less energy) than that of the exciting photon. s V* ci ™ atld lifetimes of components in a mixture. See 
, t . , , . ,. « t , # Burdick, D. S., et al., Resolution of Multicomponent Fluo- 
If we apply a 'monochromatic', short light pulse to the \t. » . \ . . f . „ . . y „ . . 
, tea . • i 4 t« . .j |i 50 rescent Mixtures by Analysis of the Excitation-Emission- 
molecules of a fluorescent material at a time to, ideally \ nu a k oo/ioom i 
represented bv Frequency Array, J. Chemomelncs, 4, 15-28,(1990). In 
" addition to requiring a non-linear minimization, a disadvan- 
tage of this method is that it requires evaluation over a large 
E(t)-E 0 -b(t-td (l) (3 dimensional) parameter space, 
where E 0 is the total energy in the pulse and E(t) is the 55 J** P reseQt ^ention utilizes a novel system that uses a 
intensity of the incident pulse as a function of time, we find chemometric analysis for measuring the fluorescence life- 
that the fluorescent intensity at a fixed wavelength has the time and lhe fluorescence spectrum of a target mixture. 
f omv ' Instead of irradiating the target sample with a single short 

pulse of light (photon counting) as other prior art systems, 
60 the present invention continuously irradiates the target mix- 

/(/Woe - *"* 0 ^ (2) mre v^i^h a ijght source whose amplitude (intensity) and 

w here phase (frequency) are modulated (varied) with time. In 

evaluating experimental phase resolved fluorescence spec- 
troscopy (PRFS) data, the variables sought are the 

/oPa-£ 0 (3) 65 conccntra ti ons> C{9 the lifetimes, t,-, and the emission spectra, 

and a has appropriate units. Thus, by measuring the decay a-fo), of the individual fluorescent contributors. The vari- 

of the fluorescent intensity with time and fitting a straight ables that we may control include the emission wavelength, 
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X, the modulation frequency, co, and the detector reference From which we see that G and S are proportional to the 
phase angle ® R . We first consider an amplitude modulated Fourier transform of Equation (11). 
excitation: In the case of noiseless data, where r fluorophores are 

present, it is necessary to make measurements at r values of 
5 X and r values of oo to obtain: 

fft-fJl+M, sin ((*)!, t>o (5) 

where w is the modulation angular frequency, M e is the g-cc-a and s-cs-a, (12) 
modulation amplitude and Eq is the average amplitude of the 

excitation. where 

At steady state, the modulated part of the fluorescent 10 
intensity may be given by 

G-{G ik \ 

Irffi^EoWA sin (aw-<l>) (6) -HS/J 

where A is the amplitude response at the fluorescent wave- 15 cc=[cos 2 <t> ( (<»/M«7,] 

length and <P is the phase lag at this wavelength and ^ sin ^ )H j 
frequency. 

We originally detected this signal (the emitted signal) by A^ftM a ik\ ( 13 ) 
correlating it with a reference signal (the driving/reference 

signal) by mixing the two signals with the mixing signal and 20 whcre the quantities defined in boldface are matrices. The 

integrating them over some multiple of a period. The cor- matrix mdices h and k run over the components, the 

related signal could therefore be expressed as: modulation frequencies, and the emission wavelengths, 

respectively. 

2jt (7) If no noise is present, we may take the indices i, j, and k 

s x (o>, = — f w M')sin(wr -<t> R )dt 25 above to run from 1 to r, the number of components. As the 

^ ^° fluorophores are assumed to be independent and as the 

= E 0 M e Aco4<t>fi -*) frequencies {coy} may be chosen arbitrarily, all the matrices 

( r defined above have true inverses. Thus, from Eqns. (12) and 

= -E 0 M e y ejcosftcosfe - ) (13) we find: 

2 U 30 

cc-^G^A (14) 

where r is the number of components present (or 
fluorophores) and <I>^ is the phase of the reference signal and 
relative to the backscattered excitation signal. We observe 35 
that Sj^oo/E^ may be written as: 

1 (8) 
S K (u>, = -E 0 M e [G(u, A)cos*ff + S(<o, A)sin4>«] ™ 

40 

Therefore, as an alternative to Eqn. (8), we may determine 
G(a),X) and S(a),X) to within a common normalization factor, Now, we note that: 
E 0 M Jl y for any fixed u) and X by making measurements at 

two different, but known, values of <b R . In the presence of « = [covKwyJsin^.-tw;) (15) 

noise, more reference phase values would be used in the 45 

measurements and fit in a least squares manner to G and S. = [co^(w;)tanv>i(w,-) 

We then find that for each value of co and X (emission 
wavelength) we have (modulo the normalization factor): 

By substitution from tan ^--oyt,-, we have: 

G(6j,A)= y £/(A)cos 2 ttM W 50 „ , vv , . 

= wccT 



S«cs-cc' u G 



SG-\ 



S(oj, A) = ^Ci(A)cosftMsiiiw(aO 



(10) 



55 



where w-diag[a) y ] and T-diagJrJ. Hence, we find 

SG-^cc-hrcc-T (17) 



60 



with tan ^<az i9 where i runs over the number of fluoro- 
phores present and assumed to be acting in an independent 
(uncorrected) fashion. 

The preferred embodiment of the present invention allows 
us to determine G and S by mixing the emitted signal from 
the mixture down to an intermediate frequency of 10 kHz, 
sample and digitize it at 50 kHz, Fourier transforming the w~ 1 s g~ u cc-cc T 

result, and picking off the signal at 10 kHz. When this is c a i / 

done the raked signal has the form: . ^ 6C V all0tt 15 10 the u forD ? of a stimda f d Rvalue/ 

65 eigenfunction equation, where the x t are the eigenvalues and 
the columns of cc are the eigenfunctions. We may also write 
I k Q)-EcAfAG sin (to //r /)-5 cos (w,^)] (li) this in the form: 



5,981,957 

12 

number of components. Furthermore, some degree of noise 
will be present in the experimental data. Therefore, it is 
useful to determine a means of solving for lifetimes and 
\ spectra when the number of frequencies and the number of 
J 5 emission wavelengths are unequal and when each is greater 
than the number of independent components present in the 

(19) mixture. 
In the following, we shall therefore assume that r com- 
ponents are present in a mixture, for which we have taken 

10 data at N frequencies and M emission wavelengths. Wc 

(20) assume that N,M^r and that N*M, in general. 
From Eqns. (12) and (16), we therefore have: 



11 



From which we have: 

r 

i-i 

and 

r 

i-i 

Since 4>X<i)/) k known once x 4 - is known, we may determine 15 
cc or cs and solve the equation: 

G=ccA or S=csA (21) 

20 

for A with a standard LU decomposition, which gives us 

A-r^ccpiw-iS (22) 

where ccpi is the pseudo inverse of cc. 25 

Another, equivalent, way of writing the eigenfunction 
equations proceeds from the observation that: 



G=ccA (28) 
(NxM)-(Nxr)-(rxM) 

and 



S-wccTA, (29) 
(NxM)<NxN)-(Nxr)-(rxr)-(rxM) 

where the number of rows and columns of each matrix is 
indicated below it. In the noise free case, all the matrices are 
of rank r. In the presence of noise, G and S will be of rank 
N or M, whichever is smaller. 

In the absence of noise, all matrices are of rank r. We 
manipulate Eqns. (28) and (29) slightly, to write: 



GA^-cc (23) 30 

and c-ccM 

and 

S-csA-wccTA (24) 35 

SO that w- l -S-cc-T-A (30) 

where I is the rxr identity matrix. 
w~ l -S-A' 1 'ce-T'G'A' l 'T (25) We may now make a singular value decomposition on the 

40 first of Eqns. (30): 

or 



G^w^SA-'^T (26) 



G-UCi'V T (31) 
(NxM)-(Nxr)-(rxr)(rxM) 

Solving this eigenequation determines x f - as the eigenval- 45 
ues and A -1 as the eigenvectors. However, as the eigenvec- where c i * the diagonal matrix of singular values of G and 
tors are normalized, A may be better determined as above if u and v are each column orthonormal, i.e.: 
information about the component concentrations is required. 

Once the spectrum and lifetime of each component is r/^iM^-v-/ (32) 

known, we may use the information to identify the fluoro- 50 

phorc. The concentration of each fluorophorc may then be From Eqns. (31) and (32), we have: 
determined by comparing its spectrum against its standard 
spectrum. Wc find: 

q e ; (A)/a ( (A)r ; (27) 



55 



with 



C,-U r GV (33) 

We define: 

Ci-V T w- l SV (34) 



60 where C 2 is not diagonal, in general, but is rxr. 

It may easily be shown, by a process that is similar to 
Gram-Schmidt orthogonalization, that while, in general, 
where a£k) is the standard spectral value of the ith com- neither matrix ^ inverlib ie, cc has a left inverse and A has 
ponent at wavelength X and t 4 . is its lifetime. a ■ ht mverse: 

The number of components in a mixture may not be 65 
known a priori, and the method requires the number of 

frequencies and wavelengths employed to at least equal the ccii coi and a ah~i (35) 
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From Eqns. (30) and (31), we may therefore write: A second method for solving Eqns. (30) is to use Eqn. (31) 



cc-UPu and A-P V >V T (36) 

with 5 

Pj^C^Ari and P^cdi UCy 



to write the pseudo inverse of G: 

Gpi-VC- x U T (43) 

so that we have: 

GGpi-UU r and GpiG-V>V T , 

We may interpret the first of Eqns. (36) in terms of its 10 where UU r is and NxN matrix and V V r is an MxM matrix, 
column vectors as the representation of the r vectors of cc in an d eac h 0 f these matrices is of rank r. 
the orthonormal basis of the r vectors of U: From Eqns. (30) and (35), we may write: 

v- (37) 
<Lj=ljU. t Pij 15 GAri-cc (44) 



i=l 



Substituting this into the first of Eqns. (30), we have: 



where c y is the jth column of cc, U,- is the ith column of U 

and P iy is the corresponding element of P^. From this 2Q Gpiw- l SAri-vv T AriT (45) 

expression and the linear independence of the c, and U.-, we , 

see that P„ is invertible. Linear independent of the c y From W and < 36 > we have: 

requires that: 



(38) 25 

2_j<*iCj = 0 only if a, = 0 for aJi j = l r. From Eqn. (32), we therefore find: 



w- l S-UPuTP v T. (46) 



From Eqn. (36), we see that Eqn. (37) implies: 



30 Substituting this into Eqn. (45), we therefore have: 



Gjxw- l 'SVV T Ari=VV r Ari'T (47) 

Eqn. (47) is the eigenvalue equation for the matrix, 
From the linear independence of the U y , we must have 35 Gpi-w -1 -S, with eigenvalues equal to the diagonal elements 

of T and eigenvectors proportional to the columns of 
(39) V V r Ari. The spectra, A, may be found as described above. 

^ Pijatj =0, for all i m l r. When determining the lifetimes, we can renormalize Eqn. 

^ (47) by using 

40 

Thus, if the condition of Eqn (37) is to be satisfied, the pscudoinv^) w l s cvec-*vec T (48) 

system of Eqns. (39) must have only the trivial solution. This 

is equivalent to the requirement that where evcc=V V r Ari for T=diag(x 1 , . . . , t r ) 

For the case when noise is present, let us assume that r 
45 components are actually present, and let us take N<M. The 
dct [P^o. experimentally determined matrices, Gx and Sx, are there- 

Therefore P„ has an inverse. A similar argument shows for f ° f ra f «* We wish to determine the NxM matrices, G 
that P is also invertible ' r r minimize some suitably defined error. 

From Eqns. (30), (33)' and (34), we have: cn We sha11 define the Euchdean norm of a malrix » M > t0 be 

50 given by: 



C 2 -Pu TP v and C x -PuIP v . (40) 

Therefore, we may write: 

C X - l C 2 ^Py l 'TPy (41) 

or 

C- l CiP v - l -P v - l T. (42) 

Eqn. (42) is the eigenvalue equation for the matrix, 
C J " 1, C 2 , with eigenvalues equal to the diagonal elements of 
T, and eigenvectors proportional to the column vectors of 
Py." 1 . The spectrum, A, may be determined from Eqns. (13) 
and (30) since cc is determined once the lifetimes are 
known. 



IIMIJ 2 = T«M T • M) = T«M • M T ) = J^lmJ 2 < 49 > 

55 

The error we define will depend on our knowledge of the 
statistics of the experimental data. We shall assume that the 
errors at each data point are uncorrelated with those at other 
data points and that the standard deviations from the sample 

60 means are the same for all data points. This is tantamount to 
the assumption that the statistics are independent of the 
modulation frequency and the emission wavelength for the 
ranges of the variables explored. Previous experience with 
experimental data tends to bear this assumption out. 

65 Given these assumptions, it is reasonable to define the 
error, x 2 , in terms of the equally weighted distance of the 
experimental points from their analytical counterparts, i.e.: 
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X 2 -||Gx-C|p + ||w- l -5r-»v- l -5|p 



16 

and Mu^Gx^QuGx+iw-^SxY'Qaiw'^Sx) 



(66) 



(50) 



where G satisfies Eqn. (31) and where, from Eqns. (40) and 
(46), we see that S satisfies: 



w- l s-uc,v T 



(5J) 



We note that M v is an NxN matrix, M^is an MxM matrix, 
and both matrices are symmetric. 

If we assume that V is known, we may determine U by 
requiring that it be chosen to maximize Eqn. (63) subject to 
the constraint imposed by Eqn. (32). We may solve for U by 
means of Lagrange multipliers: 



Therefore, from Eqns. (31), (50) and (51), we must find U, 10 Define: i\'~r\-Ti[L v iu T -u-f)] (67) 

V, C A and C, to minimize: 

where Ly is a symmetric rxr matrix formed from the 
multipliers. 

X 2 ^|G*-t/■c 1 •v r7 F+l|H'- , -Sx-UCtVT (52) We may now extremize if with U unconstrained to find: 

15 

Minimizing Eqn. (52) will yield a least squares fit for the 
data. MyU=ULy (68) 

For Eqn. (52) to be a minimum, its partial derivatives with . . . T , - , 

respect to each free variable must vanish. Applying this . S ™ lt 15 symmetric, L v may be diagonahzed by a 
requirement with respect to the variables of C lf we find: ^ transform wth an appropriate orthogonal matrix. 



C^U^Gx-V 

or, from Eqn. (31), 

G=UU T <jrVV r =Q u GrQ Y 

where 

Qu-U U 7 and Qy=VV T 
Proceeding similarly for the variables of C 2 , we find: 

C 2 -U T w- l >SrV 

or w-i-S-QvW^Sx Qv 

or mUlF-w^ S-Vx-V* 



(53) 



Therefore, we may write: 



m v v=uodo- x 



(69) 



(54) 



25 where D is an rxr diagonal matrix and 0 is an rxr orthogonal 
matrix. 
We may then write: 



30 



(55) 



35 

(56) 
(57) 

(58) 40 



MyUO-U-O-D 

Substituting Eqn. (70) into Eqn. (63) we find: 

j) = Tr{My 'U'O'O' 1 -U T ) 
= Tr[(U'OyD-(UOf) 



(70) 



or, finally, 



Eqns. (54) and (57) determine C 1 and Q once U and V are n-M^)- ( 7 °) 

determined. We now proceed to determine these variables. 

By substituting for G and w _1 S from Eqns. (54) and (57) Thus, to maximize tj we choose U to maximize D as 
into Eqn. (50), we find that minimizing that expression 45 determined by Eqn. (70). We may easily do this by observing 

that Eqn. (70) has the form of an eigenvalue equation. The 
matrix D yields r of the N possible eignevalues of the matrix 
M v . We therefore solve the eigenvalue equation for M v : 

(59) 



implies that we must maximize: 

T}M\Qu-Gx-GMMQS» l 'Sx-QX 



or 

Observing that: 

QirQu* QuQu-Qu> 
and Q v Qv-Qv, 

we may rewrite Eqn. (59) in the alternative forms: 

r\-Tr{M v 'U'Ur) 
or n-TKAft/V V 7 ) 

where My^xQ v Gx T +{w- l SxyQy {w- l 'Sx) T 



50 



MyO-OE 



(71) 



(60) 



taking the largest r eigenvalues of E for D and, without loss 
of generality, choosing U to equal the corresponding eigen- 
vectors of U. 

55 A similar procedure may be used to determine V from 
Eqn. (64) if U is known using the relationship M^/V-V-L^ 



(61) 

(62) 



FIG. 8 is a block diagram of the preferred embodiment of 
the present invention that uses the above foundation as a 
60 basis for a chemometric analysis of the emitted signal to 
determine the fluorescence spectrum and fluorescence life- 
times. The present invention first acquires the data from the 
data signal and the reference signal to determine the fre- 
( 63 ) quency and wavelength vectors 132. Next, the present 
65 invention assumes an arbitrary but fixed number of 
fluorophores, r, 134. We then determine the initial G and S 
(65) matrices 136. 



(64) 



X>2A^ (76) 
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As the frequencies in w are large, we renormalize Since we collect an amplitude and a phase value for each 

||w~\Sx|| to equal ||Gx|| so as to avoid skewing the results in emission wavelength and each modulation frequency, we 

favor of the Gx data 138. This amounts to a resealing of the have: 
units for the frequencies and the lifetimes. This resealing of 
units is compensated once the lifetimes are found, so they 5 
are expressed in seconds. 

We begin the iterative approach by determining the initial where N w is the number of frequencies and N x is the number 

approximations for U, V, g and w -1 s approximations 140, of wavelengths used. 

142, and 144. The initial iteration begins by choosing the r From Eqns. (54), (57) and (68), we see that the free 

columns of \J 1 to be the eigenvectors corresponding to the 10 parameters being fit to the experimental data are the ele- 

r largest eigenvalues of the matrix ments of U, V, Ly and h v subject to the constraints of Eqns. 

(32). For r components, we therefore find that: 

Gx-Gx r +(w- 1 -Sx)-(w- 1 -SjE) r (73) ^ ^ ^ 

or [GG^w-^SXw-'.SiT-U-UL (74) 15 Nf = r-N w +r-N x +2 -Hr+ l)-2. --r(r+ 1) 

where L=diag(>^, . . . , X r ) for ^> . . . >X r . =r-N w + r/v A 
Next, V 1 is used to form the matrix in accordance with 
Eqn. (66). We can determine W 1 from M m by solving the 

eigenfunction equation V^V Ly by choosing its r col- w ^ evaluate % 2 r for eacn r > and choose the value of r 

umns to be the eigenvectors corresponding to the r largest 20 corresponding to the smallest value of x 2 r to represent the 

eigenvalues of M w . We then form M M in accordance with number of components which best fits the data 158. 

Eqn. (65). If two or more lifetimes are the same, their components 

I s * We next initialize the decision variables 146 to produce will not be separable by this method, nor will they be 

O the solution in a reasonable time, with a reasonable accuracy, separable by the method of PREEMs. Their contributions 

]s» and a reasonable number of iterations. We begin the iterative 25 will remain combined as that of a single component with the 

!"{ solution by comparing the decision variables 148 for con- common lifetime and a spectrum given by the combined 

iJj vergence to see if further determinations for U, V, g and w _1 s spectrum appropriate to their concentrations. 

|Jl are necessary as above. In the iterative solution, we proceed FIG. 9 illustrates an embodiment of the present invention 

p to maximize r\ 7 given by Eqn. (59), and recognizing that: for the method described in FIG. 8 for a chemometric 

~ 30 analysis of the emitted signal to determine the fluorescence 

fjl 2 spectra and fluorescence lifetimes. The user interface 170 

0SnsyGx|f+l|w- 1 .Sx|| -N (72) allows the user of the present invention to interact with the 

Proceeding in this way, we form a sequence of r| values system of the present invention. The user interface allows 

;i {r| t }. We see easily that the sequence is monotonic non- the user to choose the number of components or fluoro- 

decreasing. By Eqn. (72) and the Bolzano- Weierstrass 35 phores to solve for. The data acquisition component 172 uses 

1^ theorem, the sequence has at least one limit point in the the processor 50 coupled to a data gathering component 

JgH- interval [0,N]. The monotonic nature of the sequence guar- through a data input/output path 68 to acquire the data from 

antees that there can be but one limit point. Therefore the the emitted signal to determine the frequency and wave- 

=j3 sequence must converge. length vectors. A determiner 174 determines the G and S 

i«3 The values of U and V corresponding to the limit of the 40 matrices from the frequency and wavelength vectors. A 

1^ sequence are then used to form G and w _1 S in accordance renormalizer 176 then renormalizes w" 1 S JC to G^. A deter- 

l* 8 with Eqns. (54) and (57). Eqn. (47) may subsequently be miner 178 next determines the initial U approximation, and 

solved for the lifetimes and spectra as described previously a determiner 180 determines the initial V approximation, 

160, 162, and 164. and a determiner 182 determines the initial g and w _1 s 

In the preceding, we assumed an arbitrary but fixed value, 45 approximations. An initialize r/updater 183 then initializes 

r, for the number of fluorescing components present. Fits to the decision variables. The iterative solution of the present 

the data may be found for different values of r. The question invention uses a comparator 184 that compares the decision 

then arises as to the value of r that best fits the experimental variables for convergence to see if further determinations for 

data. Two qualitative criteria may be stated: (1) the spectra U, V, g and w" J s are necessary. If the solution does not 

should look reasonable; and (2) the lifetimes should be 50 converge, the present invention uses the determiner 178 to 

positive real numbers. determine the next U approximation, the determiner 180 to 

Somewhat more quantitatively, we recognize that x 2 » determine the next V approximation, and the determiner 182 

suitably renormalized for frequency as described above, to determine the next g and w -1 s approximations. The 

represents our error for a fixed value of r. However, the value present invention then uses the initiahzer/updaler 183 to 

of x 2 depends on the number of free parameters available for 55 update the decision variables. 

the fit, which, in turn, depends on r. Instead, we attempt to When the solution converges, a determiner 186 deter- 

remove this dependence by dividing by the difference mines the reduced best fit error for the analysis of the 

between the number of experimental data points and the individual fluorescence lifetimes and fluorescence spectra, 

number of free parameters present: The present invention then uses an extractor 188 that 

60 extracts the individual fluorescence lifetimes from the target 

f (75) mixture, and an extractor 190 that extracts the individual 

= n x - JT f fluorescence spectra. 

FIG. 10 illustrates the present invention's ability to iden- 
tify and discriminate between individual overlapping spec- 
where x 2 r is the reduced error term, Ny is the number of 65 tral components in a target mixture. The present invention 

experimental data points and N f is the number of free allows for the extraction of both the individual fluorescence 

parameters. spectra and lifetimes from a target mixture of fluorophores. 



